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■^ Abstract 

i-G We study by a combination of numerical and analytical Evans function techniques 

C^ the stability of solitary wave solutions of the St. Venant equations for viscous shallow- 

G water flow down an incline, and related models. Our main result is to exhibit examples 

1—^ of metastable solitary waves for the St. Venant equations, with stable point spectrum 

indicating coherence of the wave profile but unstable essential spectrum indicating os- 

►^ dilatory convective instabilities shed in its wake. We propose a mechanism based on 

P^ "dynamic spectrum" of the wave profile, by which a wave train of solitary pulses can 

\^ stabilize each other by de-amplification of convective instabilities as they pass through 

OQ successive waves. We present numerical time evolution studies supporting these con- 

'^ elusions, which bear also on the possibility of stable periodic solutions close to the 

t^^ homoclinic. For the closely related viscous Jin-Xin model, by contrast, for which the 

^P essential spectrum is stable, we show using the stability index of Gardncr-Zumbrun 

that solitary wave pulses are always exponentially unstable, possessing point spectra 

with positive real part. 
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1 Introduction 

Roll waves are a well-known phenomenon occurring in shallow water flow down an inclined 
ramp, generated by competition between gravitational force and friction along the bottom. 
Such patterns have been used to model phenomena in several areas of the engineering liter- 
ature, including landslides, river and spillway flow, and the topography of sand dunes and 
sea beds and their stability properties have been much studied numerically, experimentally, 
and by formal asymptotics; see |BMj and references therein. Mathematically, these can be 
modeled as traveling- wave solutions of the St. Venant equations for shallow water flow; see 
[D1 IN1[ IN2j for discussions of existence in the inviscid and viscous case. They may take the 
form either of periodic wave-trains or, in the long-wavelength limit of such wavetrains, of 
solitary pulse-type waves. 

Stability of periodic roll waves has been considered recently in |NHIN2[IJZNj . Here, we 
consider stability of the limiting solitary wave solutions. This appears to be of interest not 
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only for its relation to stability of nearby periodic waves (see |GZl lOZlj and discussion in 
|JZNj ). but also in its own right. For, persistent solitary waves are a characteristic feature 
of shallow water flow, but have been more typically modeled by dispersive equations such 
as Boussinesq or KdV. 

Concerning solitary waves of general second-order hyperbolic-parabolic conservation or 
balance laws such as those examined in this paper, all results up to now [AMPZl] IGZ|, IZ2| 
have indicated the such solutions exhibit unstable point spectrum. Thus, it is not at all 
clear that an example with stable point spectrum should exist. Remarkably, however, for 
the equations considered in |N2t IJZNj we obtain examples of solitary wave solutions that 
are metastahle in the sense that they have stable point spectrum, but unstable essential 
spectrum of a type corresponding to convective instability. That is, the perturbed solitary 
wave propagates relatively undisturbed down the ramp, while shedding oscillatory instabil- 
ities in its wake, so long as the solution remains bounded by an arbitrary constant. The 
shed instabilities appear to grow exponentially forming the typical time-exponential os- 
cillatory Gaussian wave packets associated with essential instabilities; see the stationary 
phase analysis in |OZlj . Presumably these would ultimately blow up according to a Ricatti 
equation in the standard way; however, for sufficiently small perturbation, this would be 
postponed to arbitrarily long time. For the time length considered in our numerics, we can 
not distinguish between the linear and nonlinear growth. 

We confirm this behavior numerically by time evolution and Evans function analysis. 
We derive also a general stability index similarly as in \G7^, IZ21 IGol IZ4] counting the parity 
of the number of unstable eigenvalues, hence giving rigorous geometric necessary conditions 
for stability in terms of the dynamics of the associated traveling-wave ODE; see Section 



3.3, Using the stability index, we show that homoclinics of the closely related viscous 



Jin-Xin model are always unstable] see Section 3.4 We point out also an interesting and 
apparently so far unremarked connection between the stability index, a Melnikov integral for 
the homoclinic orbit with respect to variation in the wave speed, and geometry of bifurcating 
limit cycles (periodics) , which leads to a simple and easily evaluable rule of thumb connecting 
stability of homoclinic orbits as traveling wave solutions of the associated PDF to stability 



of an enclosed equilibrium as a solution of the traveling-wave ODF; see Remark 3.13 

The examples found in this paper are notable as the first examples of a solitary-wave 
solution of a second-order hyperbolic-parabolic conservation or balance law for which the 
point spectrum is stable. This raises the very interesting question whether an example could 
exist with both stable point and essential spectrum, in which case linearized and nonlinear 
orbital stability (in standard, not metastable sense) would follow in standard fashion by the 
techniques of |MaZl[ IMaZ4[ [Z2l iJZl iJZNl ILETZl ITZT] . In Section [5] we present examples 
for various modifications of the turbulent friction parameters (different from the ones in 
|N2t IJZN| ) that have stable point spectrum and "almost-stable" essential spectrum, or 
stable essential spectrum and "almost-stable" point spectrum, suggesting that one could 
perhaps find a spectrally stable example somewhere between. However, we have up to now 
not been able to find one for the class of models considered here. Whether this is just by 
chance, or whether the conditions of stable point spectrum and stable essential spectrum 



1 INTRODUCTION 4 

are somehow mutually exclusivqjis an extremely interesting open question, especially given 
the fundamental interest of solitary waves in both theory and applications. 

Finally, we note that similar metastable phenomena have been observed by Pego, Schnei- 
der, and Uecker |PSU] for the related fourth-order diffusive Kuramoto-Sivashinsky model 

Uf + o^u + o^u -\ — = U, 

an alternative model for thin film flow down a ramp. They describe asymptotic behavior of 
solutions of this model as dominated by trains of solitary pulses, going on to state: "Such 
dynamics of surface waves are typical of observations in the inclined film problem |CD] , both 
experimentally and in numerical simulations of the free-boundary Navier-Stokes problem 
describing this system." 

That is, both our results and the results of |PSUj seem to illustrate a larger and some- 
what surprising phenomenon deriving from physical inclined thin-film fiow of asymptotic 
behavior dominated by trains of pulse solutions which are themselves unstable. We discuss 
this interesting issue, and the relation to stability of periodic waves, in Section [6j suggest- 
ing a heuristic mechanism by which convectively unstable solitary waves, when placed in 
a closely spaced array, can stabilize each other by de-amplification of convected signals as 
they cross successive solitary wave profiles. We quantify this by the concept of "dynamic 
spectrum" of a solitary wave, defined as the spectrum of an associated wave train obtained 
by periodic extension, appropriately defined, of a solitary wave pulse. In some sense, this 
notion captures the essential spectrum of the non-constant portion of the profile. This is 
in contrast to usual notion of essential spectrum which is governed by the (often unstable) 
constant limiting states. For the waves studied here, the dynamic spectrum is stable, sug- 
gesting strongly that long- wave periodic trains are stable as well. This conjecture has since 
been verified in |B.TNRZl[lB.TNRZ2j . 

Ackno'wledgement: Thanks to Bjorn Sandstede for pointing out the reference |PSU] . 
and to Bernard Deconink for his generous help in guiding us in the use of the SpectrUW 
package developed by him and collaborators. K.Z. thanks Bjorn Sandstede and Thierry 
Gallay for interesting conversations regarding stabilization of unstable arrays, and thanks 
the Universities of Paris 7 and 13 for their warm hospitality during a visit in which this 
work was partly carried out. The numerical Evans function computations performed in this 
paper were carried out using the STABLAB package developed by Jeffrey Humpherys with 
help of the first and last authors. 



^For a similar dichotomy in a related, periodic context, see |0Z1| . 
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2 The St. Venant equations with viscosity 

2.1 Equations and setup 

The 1-d viscous St. Venant equations approximating shallow water flow on an inclined 
ramp are 

ht + (hu)^ = 0, 
{hu)t + {h^/2F + /in^)^ = h- u\u\''-^/h' + u{hu^)^, 

where l<r<2, 0<s<2, and where h represents height of the fluid, u the velocity 
average with respect to height, F is the Froude number, which here is the square of the ratio 
between speed of the fluid and speed of gravity waves, u = Re~^ is a positive nondimensional 
viscosity equal to the inverse of the Reynolds number, the term u\u\^~^ /h'^ models turbulent 
friction along the bottom, and the coordinate x measures longitudinal distance along the 
ramp. Furthermore, the choice of the viscosity term v{hux)x is motivated by the formal 
derivations from the Navier Stokes equations with free surfaces. Typical choices for r, s 
are r = 1 or 2 and s = 0, 1, or 2; see [BMl INH IN2j and references therein. The choice 
considered in [Nil |N2l llZN] is r = 2, s = 0. 



Following [JZNj . we restrict to positive velocities u > and consider (2.1 ) in Lagrangian 
rd 

(2.2) 



coordinates, in which case (2.1) appears as 

Tt-Ux = 



ut + ((2F)- V-2), = 1 - t'+^u'- + v(t- 



Ux 



where r := h^^ and x now denotes a Lagrangian marker rather than physical location. 
Note that working in Lagrangian coordinates was crucial in developing the large-amplitude 
damping estimates necessary for the nonlinear stability analysis of periodic roll waves in 

[JZNJ. 

Remark 2.1. Since we study waves propagating down a ramp, there is no real restriction 
in taking n > 0. However, we must remember to discard any anomalous solutions which 
may arise for which negative values of u appear. 



We now consider the existence of solitary traveling wave solutions of (2.2). Denoting 



U := {t,u), consider a traveling-wave solution U = U{x — ct), of (2.2), which is seen to 
satisfy the traveling-wave ODE 

-cf' -u' = 0, 
^^■^^ -cu' + ((2F)-ir-2)' = 1 - f'+^u' + i^{f-^u'y. 



Integrating the first equation of (2.3) and solving for u = u{t) ■= q — ct, where q is the 



resulting constant of integration, we obtain a second-order scalar profile equation in r alone: 

(2.4) c^f' + {{2F)-^f-^y = 1 - f^+i(g - cff - cvif-'^f') . 
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Note that nontrivial traveling-wave solutions of speed c = do not exist in Lagrangian 
coordinates, as this would imply u = q, and (2.4) would reduce to a scalar first-order 
equation 



(2.5) 



Fr\f'+'q' - 1 



which has no nontrivial solutions with r > 0. For non-zero values of c, however, it is seen 



numerically that (2.4) admits homoclinic orbits corresponding to solitary wave solutions of 



(2.2). See also |Nlt IN2j . Furthermore, generically these are seem to emerge from a hopf 



bifurcation from the equilibrium state, which generates a family of periodic orbits which 
terminate into the bounding homoclinic. 

2.2 The subcharacteristic condition and Hopf bifurcation 



To begin, we consider the stability of the equilibrium solutions (2.2). First note that (2.2) 
is of 2 X 2 viscous relaxation type and can be written in the form 

(2.6) Ut + H{U)^-v{B{U)U^),-. 

where g{U) = 1 — t^^'^u^' and 





9{U) 







H{U)=[^\, B{u) = 

\2Ft2/ 

At equilibrium values u = t~'''+^'''" > 0, the inviscid version of (2.6), obtained by setting 

of dH, and 



in (2.6), has hyperbolic characteristics equal to the eigenvalues it- 



-3/2 



equilibrium characteristic 



's+l^ 



-{r+s+l)/r 



equal to df^{T), where /*(t) := g{T,u^{T)) and 



u^,{t) := T~^^~^^>'^ is defined by g{T,u^{T)) = 0. The subcharacteristic condition, i.e., the 
condition that the equilibrium characteristic speed lie between the hyperbolic characteristic 
speeds, is therefore 



(2.7) 



dUr) 



s + l\ u^{t) ( s + 1 



-{r+s+l)/r 



-3/2 



< 



F 



Cs{t). 



For the common case r = 2, s = considered in |N2| IJZNj . this reduces to F < 4. For 
2x2 relaxation systems such as the above, the subcharacteristic condition is exactly the 
condition that constant solutions be linearly stable, as may be readily verified by computing 
the dispersion relation using the Fourier transform. 



Remark 2.2. Let q and c be fixed. Then the equilibria of the traveling-wave ODE (2.4) are 
given by zeros of g^f^r) = g{T, q — ct) (corresponding to u^:{t) = q — ct). Notice the relation 
dg*{j) = gr{T) — cguir) = guiT){df^:{T) — a). Since g^ < 0, by the relaxation structure of 
( |2.6[ ), the key quantity df^^ — c has a sign opposite to the one of dg^,. Therefore, provided 
there is no multiple root of g^,, the sign of dfi,{To) — c alternates among the equilibria tq of 



the traveling-wave ODE (2.4). 
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For the full system (2.6) with viscosity z^ > 0, we now show some Fourier computations 



relating the subcharacteristic condition with both stability and Hopf bifurcation. Let speed 



c be fixed. Then rewriting (2.6) in a moving frame {x — ct, t) and linearizing about a constant 
solution (r, n) = (To,no), with tq > and uq = w*(to) > 0, we obtain the linear system 
Ut + AoUx = BQU + CoUr,^, wherfj3 



(2i 



and c. 



^n 



-c 

^2 



Bn 





-(s + 1)/to 





-r/uo 



Co 







UT, 







-3/2 



Cs[To) = "r^ denotes the positive hyperbolic characteristic of the inviscid problem 
(see above). The corresponding dispersion relation between eigenvalue A and frequency k 
is readily seen to be given by the polynomial equation 



(2.^ 



= det(AI + Aoik - Bq + k^Co) 

vk^' 



\' + 



Uq 



2ick + 



X + ik 



s + 1 



cr 



Uq 



+ ik{c^ 



cvk^ 



Inspecting stability, we Taylor expand the dispersion relation (2.9) about an eigenvalue A 



with 3f?(A) = 0. Note that if A is such an eigenvalue corresponding to a non-zero frequency 
k, we may assume A = by replacing c with c — 9(A)/A;. Let us then focus on A; = or 
A = and /c / 0. 

Eignevalues corresponding to /c = are and —r/uQ. Solving the dispersion equation 



(2.9) about (A,/c) = (0,0) with A(/c) and differentiating, one finds 

s + r 



and 



-A"(0) 



no 
r 



A'(0) 



(iA'(0) + c)' 



no 
To 



liO 
r 



s + 1 
r 



no 
To 



so that failure of the subcharacteristic condition implies spectral instability. Furthermore, 
it follows that any solitary wave solution of (2.2) for which the limiting value lim^^-s^ioo ''"(x) 



violates the subcharacteristic condition (2.7) will have unstable essential spectrum 



(2.10) 



We now look for eigenvalues A = for /c 7^ by solving 

s + 1 cr .. , o o. cvk'^ 







+ ik{c^ 



To no 
Such a k exists if and only if c = Cg and {^-^) ^ > c^, and then the solutions are ife/^ with 



k 



H 






s + 1 



no 
To 



Recall that r^ 



Wo 
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Now, solving the dispersion equation (2.9) about (A, A;) = (0,ibA;//) with X{k) and differen- 
tiating, one finds 



(2.11) 



X\±k 



2c.u 



H 



ik^ 



'^0 ^/^o + ^kjj/r^ =F 2iCskH 



yielding again instability. See Figure [T] for a typical depiction of the triple degeneracy of 
the essential spectrum at A = for an unstable constant state when c = Cg- 



0.08 - 

0.06 - 

0.04 - 

0.02 - 
- 
-0.02 - 
-0.04 - 
-0.06 - 
-0.08 - 

-0.1 L^ 

-0.1 



.••-) 



T= 3.927 

c= 0.57052639 



Figure 1: The essential spectrum for the unstable equilibrium solution {tq,uq) of (2.4) at 
a Hopf bifurcation c = c^ is shown, where here we take uq = 0.96, q = uq + c/u^, v = 0.1, 
r = 2, s = 4, F = 6, and c = Cg = 0.57052639. Notice the triple degeneracy of the A = 
eigenvalue in this case. 



Actually equation (2.10) is the linearization of the profile equation (2.4) (with q = 
uq + CTo) about To, expressed in Fourier variables, so that previous pragraph translates to 



the Hopf bifurcation conditions for the generalized St. Venant equations (2.2) 



(2.12) 



Csiro) 



.-3/2 



and (i/*(ro) > Cs(ro) 



Moreover k}{ is the limiting frequency associated to the Hopf bifurcation and relation (2.11) 



will yield instability for close-by periodic travelling-waves generated through the Hopf bi- 
furcation. See |JZNj for an alternate derivation of these bifurcation conditions. 



Remark 2.3. Experiments of [NU IBJNRZH IBJNRZ2J in case r = 2, s = indicate that, 
when F > 4, there is a smooth family of periodic waves parametrized by period, which 
increase in amplitude as period is increased, finally approaching a limiting homoclinic orbit 
as period goes to infinity; see |N2l IBJNRZTl IBJNRZ2] . Likewise, in the study of a related 
artificial viscosity model with r = 2, s = jHCj . periodic and homoclinic solutions were 
found in the regime of unstable constant solutions, F > A. More precisely, in [HCj this 
very configuration and generation of the periodics and homoclinics were found analytically 
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near where the nonhnear center and saddle point coalesce. The authors were then able to 
analytically follow these configurations to obtain a global description of phase space. 

Remark 2.4. Let us further examine the linearized profile equation. For this equation, if 
c((i/*(ro) — c) < then the equilibrium tq is a saddle point ; whereas if c(d/^,(ro) — c) > 
then if c{c^ — Cs(ro)^) < the equilibrium tq is a repellor, either of spiral or real type, 
and if c{c^ — Cs{tq)'^) > the equilibrium tq is an attractor, either of spiral or real type. 
In particular, by the alternating property, the repellor case describes behavior at a single 
equilibrium ri trapped within a homoclinic (asymptotic to tq) if and only if there stands 
€{0^ — Cs(ri)^) < 0. For the example wave considered in Section [s] we have r = 2, s = 0, 
c = .78, F = 6,q = 1.78, tq = 1 and Cs(ro) = 1/V2F < .289 < .78 at the saddle equilibrium, 
but, at the enclosed equilibrium n < 1/2, Cs(ri) = t^ jylF > .816 > .78. 

3 Analytical stability framework 

3.1 Linearized eigenvalue equations and convective instability 

We now begin our study of the stability of a given homoclinic solution U{x, t) = U{x — ct) 
satisfying liinx^±oo = Uq to localized perturbations in L^(]R;R^). To this end, notice that 
changing to co-moving coordinates, linearizing about U = {f,u), and taking the Laplace 
transform in time we obtain the linearized eigenvalue equations 



-2„,/\/ 



Xt - ct' -u' = 0, 
Xu-cu' - (f-3(F-i - 2uUx)Ty = -{s + l)f'uV - rf'+^u'-\ + u{f ~u , 

As we are considering localized perturbations, we consider the L^ spectrum of the eigenvalue 



equations (3.1 ). Further, we say the underlying solitary wave U is spectrally stable provided 
that the L spectrum of (3.1 ) lies in the stable half plane {3f?A < 0} and is spectrally unstable 
otherwise. 

In order to analyze the spectral stability of the underlying solution U, we consider the 



point and essential spectrum of the linearized equations (3.1 ) separately. Here, we define the 
point spectrum in the generalized sense of |PW2j (including also resonant poles) as zeros of 
an associated Evans function (description just below), and stability of the point spectrum 
as nonexistence of such zeros with nonnegative real part, other than a single translational 
eigenvalue always present at frequency A = 0. Equivalently [Satj . this notion of generalized 
point spectrum can be described in terms of usual eigenvalues of the linearized operator with 
respect to an appropriate weighted norm; see |TZ2] for further discussion. Furthermore, 
we define the essential spectrum in the usual complementary sense: that is, as the relative 
complement of the point spectrum (in the generalized sense described above) in the L^ 



spectrum of the linearized equations (3.1). 



Concerning the structure of the essential spectrum, we recall a classical theorem of 



Henry |He| IGZ| which states that the essential spectrum of the linearized operator (3.1) 



about the wave both includes and is bounded to the right in the complex plane by the 
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right envelope of the union of the essential spectra of the linearized operators about the 
constant solutions corresponding to the left and right end states lini2,_s>itoo U{x), which in 
our case agree with the common value C/q- This in turn may be determined by a Fourier 



transform computation as in Section 2.2 Consulting our earlier computations, we see that 
in the common case r = 2, s = 0, constant solutions are unstable in the region of existence of 
solitary waves, and hence solitary waves in this case always have unstable essential spectrum 



(see Remark 2.3), and hence are linearly unstable with respect to standard (e.g., L^ or H'^) 
norms. 

As discussed in |PW2| [Sat] . however, such instabilities are often of "convective" nature, 
meaning that growing perturbations are simultaneously swept away from, or "radiated" 
from the solitary wave profile, which itself remains, at least to linear order, intact. This 
corresponds to the situation that generalized point spectrum, governing near-field behavior, 
is stable, while essential spectrum, governing far-field behavior is unstable. The phenomenon 
of convective instability can be captured at a linearized level by the introduction of an 
appropriate weighted norm jSatl IPWlj . with respect to which the essential spectrum is 
shifted into the negative half-plane and waves are seen (by the Hille-Yosida Theorem) to 
be linearly time-exponentially stable. Alternatively, see the pointwise description obtained 
by stationary phase estimates in |OZlj of the far-field behavior as of an oscillatory time- 
exponentially growing Gaussian wavepacket convected with respect to the wave. (See also 
|AMPZ2j .) 

Next, we continue our discussion by developing the necessary tools to analyze the point 



spectrum of (3.1 ). Notice that, as stated in the introduction, spectral stability of both point 
and essential spectrum, should it occur, can be shown to imply linearized and nonlinear 
stability in the standard time-asymptotic sense by the techniques of |MaZll IMaZ31 IMaZ4| 
ILRTZ[ TTZU IJZN] . However, as also noted in the introduction, we have so far not found a 
case in which these conditions coexist. 

3.2 Construction of the Evans function 

In order to describe the point spectrum (defined above in the generalized sense), we outline 
here the construction of the Evans function corresponding to the linearized eigenvalue equa- 



tions (3.1). First, notice that setting W := {t,u,t u') , and recalling that u = {q — ct 



we may write (3.1) as a first-order system of the form W = A{x, X)W, which is suitable for 



an Evans function analysis, where 

-fyc 

(3.2) A:= \ r2 

—cf^+af^ /c 




(3.3) a := f-^iF'^ + 2cuf^), a^ = -3f-^f^{F-^ + 2cuf^) + 2f-^ciJf^^. 

As with the essential spectrum, the point spectrum can essentially be described by analyzing 



the behavior of the system (3.2) evaluated at the constant end state Uq of the underling 
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solitary wave. To begin, we make the following definition which extends the notion of 



hyperbolicity of the limiting first order system (3.2). 



Definition 3.1. Following I A GJ\ \U^ . for a generalized eigenvalue equation W' = vl(x, A)T^, 
with 

lim A{x,\) = Aq{\), 

we define the region of consistent splitting to he the connected component of real plus infinity 
of the set of A such that Aq has no center suhspace, that is, the stable and unstable subspaces 
S{Aq) andU{AQ) of A± are of constant dimensions summing to the dimension of the entire 
space and maintain a spectral gap with respect to zero. 

We say that an eigenvalue problem satisfies consistent splitting if the region of consistent 
splitting includes the entire punctured unstable half-plane {A 1 3f?A > 0}\{0} of interest in the 
study of stability|^ The notion of consistent splitting contains all the necessary information 
in order to define an Evans function. Furthermore, the notion of consistent splitting has an 
immediate consequence concerning the stability of the limiting end states. 

Lemma 3.2 ( [ AGJl KjZ\ IMaZ3] ). Consistent splitting is equivalent to spectral stability of 
the constant solution U = Uq, defined as nonexistence of spectra of the constant- coefficient 
linearized operator about Uq on the punctured unstable half-plane { A | KA > 0} \ {0}. 

Proof. This follows from the standard observation |GZ[ IZHl IMaZ31 IZlj that the set of 
A for which ^o(A) possesses a pure imaginary eigenvalue fi = ik is given by the linear 
dispersion relations A = Xj{k) for the associated constant-coefficient equation Ut = LU := 
BUxx - AUx + CU, defined by the roots of det(AI - k'^B - ikA + C) = 0. D 



In the common case r = 2, s = 0, the Fourier transform computations of Section 2.2 
show that constant solutions are unstable in the regime of existence of homoclinic orbitq^ 
and thus consistent splitting must fail in this case. Thus, while the notion of consistent 
splitting is sufficient for the construction of an Evans function, it is insufficient for our 
purposes. For our applications then, we need a generalization of the notion of consistent 
splitting which still incorporates the essential features necessary for the construction of an 
Evans function. Such a generalization was observed in (PW2j . which we now recall. 

Definition 3.3. We define the extended region of consistent splitting to be the largest 
superset of the region of consistent splitting on which the total eigenspaces associated with 
the stable (resp. unstable) subspaces of Aq maintain a spectral gap with respect to each 
other, that is, for which the eigen-extensions S and lA of the stable and unstable subspaces 
S{Aq) and ZY(Ao) defined for X near real plus infinity satisfy the condition that the maximum 
real part of eigenvalues of S is strictly smaller than the minimum real part of of eigenvalues 
ofU. 



^Recall, [GZI IMaZ3l IJZN] . that the special point A = may be adjoined by analytic extension. 
*At least, for those generated through a Hopf bifurcation. Recall, however, that all homoclinic orbits 
studied here are obtained in this way. 
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We say that an eigenvalue problem satisfies extended consistent splitting if the re- 
gion of extended consistent splitting includes the entire punctured unstable half-plane 



{A|KA > 0} \ {0}. Notice that the Fourier transform computations of Section 2.2 show 
that extended consistent splitting holds in a neighborhood of the origiqj in the common 
case r = 2, s = while, as mentioned above, consistent splitting does not. Next, under the 
assumption of extended consistent splitting we show that it is possible to relate the large x 



behavior of the first order system W = A{x, X)W, with A given as in (3.2), to the behavior 



of limiting system W' = Aq(X)W, where the coefficient matrix ^o is defined as above. 

Lemma 3.4 ( |GZ| \ZJH\ IZ1| ). Let A{x,X) be analytic in X as a function into C^{x), expo- 
nentially convergent to Aq as x ^ ±00, and satisfying extended consistent splitting. Then, 
there exist globally defined bases {W^ , . . . jWjT} and {W^-^, . . . , Wj^}, analytic in X on 
{3f?A > 0} \ {0}, spanning the manifold of solutions of W = A{x,X)W approaching expo- 
nentially in angle as x ^ —00 (resp. x — t- +00 j to ZY(^o) and 5(^40), where k = diioaU{Ao) 
and n = dim U; more precisely, 

(3.4) W^j^ ~ e^^o^l/^ as x ^ ±00 

for analytically chosen bases V- ofL({AQ) andS^Ao). For the eigenvalue equations consid- 
ered here, these bases (both W- and V- ) extend analytically to 3f?A < —r] < 0. 

Proof. The first, general result follows by the conjugation lemma of |MeZj (see description, 
|Zlj ) asserting that exponentially convergent variable-coefficient ODE on a half- line may be 
converted by an exponentially trivial coordinate transformation to constant-coefficient ODE 
Z' = AqZ, together with a lemma of Kato [Kj IZ6] asserting existence of globally analytic 
bases {V- } for the range of an analytic projection, the eigenprojections being analytic due 
to spectral gap. Extension to A = (and thus KA < —rj < for some rj) then follows by a 
matrix perturbation expansion at A = verifying analytic extension of the projectors. See 
\GZ\ IMaZ3[ [Zlj for the now-standard details of this argument. D 



By Lemma 3.4 then, solutions of the full linearized system W' = A(x, X)W behave 
asymptotically like the solutions of the limiting system W' = Ao{X)W. Thus, if A is in 
the point spectrum (in the usual L^ sense) of (3.1) the corresponding vector solution of 



full linearized system should decay as x — )• itcxD corresponding to a connection between the 
stable and unstable subspaces of Aq. For the linearized St. Venant eigenvalue equations 



considered here, computations like those of Section 2.2 show that on the region of consistent 
splitting the limiting coefficient matrix ^o('^) has a one-dimensional stable subspacerl and 



hence Lemma 3.4 applies with k = 2 and n = 3; as with extended consistent splitting, this 



is something that can be numerically verified away from A = 0, 00. It follows that A belongs 



^Away from the origin, extended consistent splitting can be checked numerically in the course of per- 
forming the Evans function computations outlined in Section [S] In particular, the numerical code for the 
Evans computations is written in such a way that failure of extended consistent splitting signals an error. 

^This dimension can be found near the origin or near real plus infinity and is, by definition, constant in 
the region of consistent splitting. 
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to the point spectrum of (3.1 ) provided the spaces spanjH^]^ , W2 } and span{VF3 } intersect 



non-trivially, which leads us to our definition of the Evans function. 

Definition 3.5. For the linearized St. Venant eigenvalue equations, we define the Evans 
function as 

(3.5) ^(,)^^det(I4^r,^.-,^3^)Uo 



det{V,-,Vi,K 



3 



where V- are the limiting directions of W- . Furthermore, we say A belongs to the point 



spectrum of (3.1) provided D{X) = 0. 



± 



Remark 3.6. Note that definition (3.5 ) is independent of the choice of basis elements V- 
Wj , so long as (3.4) holds. Thus, it is not necessary in practice to choose the V, analyti 



cally, or even continuously. This way of normalizing the Evans function for homoclinics by 
considering the limiting vectors V_- seems to not only be new to the literature, but also 
quite convenient from a computational standpoint; see calculations below. 



When consistent splitting holds, the basis elements W- span the manifolds of solutions 
decaying as a; — )• +00 (resp. x — )■ —00), and so vanishing of D is associated with existence 
of an exponentially decaying eigenfunction and roots A correspond to standard eigenvalues. 
When extended consistent splitting holds but consistent splitting does not, then some of 
the basis elements may be nondecaying or even growing at infinity, and so roots correspond 
rather to "resonant poles". See |PW2j for further discussion. Thus, as mentioned before, our 
definition of point spectrum is in a generalized sense including all roots (including resonant 



poles) of the Evans function as defined in Definition 3.5 However, notice that if one has 
good essential spectrum and bad point (Evans) spectrum, then the unstable Evans spectra 
lies in the region of consistent splitting (excepting the degenerate case A = 0) and hence 
corresponds to a spectral instability to localized perturbations. In particular, in this case is 
is possible to conclude nonlinear instability of the solitary wave by arguments like that of 
|Z81IZ9] . On the other hand, if one has both bad essential and bad point (Evans) spectrum, 
then the instability depends on which component of the spectra has larger real part; if 
the unstable point spectrum is farther out, then again this corresponds to an unstable 
localized eigenvalue, while if the essential spectrum is farther out, then we are essentially in 
the convective instability case briefly discussed in the introduction. The latter convective 
instability will be discussed further in Section [6} 

3.3 The stability index 

Next, we use the Evans function to seek conditions which imply the existence of unstable 
point (Evans) spectrum. First, however, we need the following lemma which guarantees that 
our construction of the Evans function in the previous section is valid along the positive 
real axis. 
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Lemma 3.7. Let tq be an equilibrium solution of the profile ODE (2.4). Then, the open 
positive real axis is contained in the set of consistent splitting for the linearized St. Venant 
equations (3.1) if and only if cdf^{To) < c^ or c^ — c^ < 0. In any case, there exists a 



neighborhood of the origin which lies in the region of extended consistent splitting. 



Remark 3.8. Notice by Remark 2.4, the condition c{df^,{TQ) — c) < is equivalent to the 
equilibrium tq being a saddle point of the profile ODE (2.4). It follows that, given any 



homoclinic solution of (2.4), the set (0, oo) is contained in the set of consistent splitting 



and hence any positive root of the Evans function corresponds to a spectral instability of the 
underlying wave to localized perturbations. 

Proof. Recall that the rightmost boundary of the essential spectrum is precisely the bound- 
ary of the region of consistent splitting. Thus, to establish the first claim it is sufficient 
to prove that the essential spectrum does not intersect the positive real axis if and only if 
cd/*(To) < c^ or c^ 



c2<0. 



This in turn follows by the Fourier analysis of Section 2.2 From (2.9) we see that that 
the essential spectrum A 



X{k) intersects the real axis if and only if 
's + 1 cr cvk"^' 



-2ckX + k 



When k ^ such roots must satisfy 
(3.6) A 



1 
2^ 



To 



s + 1 



To 



cr 

Uq 



cr 
Mo 



0. 



cv 



k^ 



and hence a sufficient condition that A < is that 



s + 1 r 



Rearranging using the fact that uq = t. 



cro 

-(s+l)/r 




UQ 



we find that this last condition is equivalent 



to ^d/=K(To) < 1, which is equivalent to c{df^,{T(j) — c) < 0. Therefore, under this sign 



condition, we conclude that the only real roots A of (2.9) that can occur when k j^ must 
be negative. But, when k 



we have the explicit roots A 
non-negative and hence yields the desired contradiction. 



and A 



Mo 



which are also 



Likewise, for k ^ 0, substituting (3.6) into the real part of the dispersion relation 2.9 
yields the relation 

A /s + 1 r uk"^ 

2 



+ 



+ 



CTq 



Uq 



kHc^ 



which is a contradiction if A > and c{df^{TQ) — c) > and c^ — c^ < 0. Thus, we also see 
that the open positive real axis lies in the region of consistent splitting whenever c^ — c^ < 0. 
If, on the other hand, c{df^{To) — c) > and c^ — c^ > 0, then, necessarily \df^\ > 
\cs\, hence the subcharacteristic condition is violated and there exist essential spectra with 
strictly real part. Fixing tq, hence df^, and Cg, now vary the wave speed c. For c = c^, we have 



3 ANALYTICAL STABILITY FRAMEWORK 15 

as already noted the Hopf configuration shown in Figure [l] of a triple degeneracy with three 
distinct roots k = 0, it/c* at A = 0, and (by the sufficient condition already established), the 
positive real axis is contained in the region of consistent splitting. Recalling that changes in 
the wave speed c by an amount Ac corresponds to a translation in the roots of the dispersion 
relation along the imaginary axis by an amount proportional to i{Ac)k, we find that for 
c^ ^ c^, the triple singularity breaks up into three distinct crossings of the imaginary axis, 
one at A = and two at ±ikit:{c — Cg). For c^ < c^, we have already shown that the real 
axis remains in the region of consistent splitting. See Figure me), which shows an example 
of the essential spectrum in a case when c > d/*(To) > c^ > 0. When c((i/*(ro) — c) > 
and cP > c^, therefore, the roots move in the opposite direction along the real axis, and the 
essential spectrum must intersect the positive real axis as in Figure ^ (Alternatively, we 
may see this by expansion of the dispersion relation in the vicinity of A = and c = c^.) 
This verifies the first claim. 



Finally, notice that expanding the dispersion relation (2.9) with the substitution /j, = ik, 
it follows that the matrix ^(x;0) = ^o(O) from the Evans system about the equilibrium 
solution (tojUq) with A = possess a one-dimensional center subspace. As the spectral gap 
between the stable and unstable subspaces of Ao{0) persist under small perturbation, it 
follows that the corresponding matrices Aq{X) satisfy extended consistent splitting for |A| 
sufficiently small, which verifies the final claim. D 



0- 

-0.02 - 
-0.04 - 
-0.06 - 
-0.08 — 



0.02 0.03 



Figure 2: This figure depicts the essential spectrum for the constant state tq ~ 1.08507 
for the common case (r, s) = (2,0) when c = 0.53, q ~ 1.53509, i' = 0.1, and F = 6. In 
particular, it is easily seen that c{df^:{To) — c) > and c^ > c^ ^ (0.361)^ in this case. Notice 
that the essential spectrum intersects the positive real axis 



Following |GZ| then, we can analyze the Evans function on a set of the form (— e, oo) 
for some e > 0. In particular, by comparing high and low (real) frequency behavior of the 
Evans function, may compute a mod-two stability index giving partial stability information 
in terms of geometric properties of the flow of the traveling- wave ODE: specifically, the sign 
of a Melnikov derivative with respect to wave speed. To this end, we begin by analyzing 
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the high (real) frequency asymptotics of the Evans function. 

Lemma 3.9. sgnD(A) — )• 1 as A — )• +oo for X real. 

Proof. This may be checked by homotopy to an easy case, followed by direct computation. 
Or, we may observe that asymptotic behavior is determined by the principal part 



A 


n 


n 


c 












f' 


* 


A 






of A{x,X), hence by the Tracking Lemma of [GZl EHl IMaZ3j . det(Tyf , 1^2"' ^3*')^=o is 
asymptotic to a positive multiple of det{V^ ,V^ ,Vr^) . See the computations of Section U^ 
for details. D 

In order to derive an instability index, we compare the above high frequency behavior 



with the low-frequency behavior near A = 0. First, notice that differentiating (2.3) with 



respect to x implies that Ux satisfies the linearized eigenvalue equations (3.1) with A = 



and hence -D(O) = 0. It follows negativity of D'{0) is sufficient for the existence of a real, 



positive, unstable element of the point (Evans) spectrum of (3.1). Notice, however, that 
-D'(O) < always implies either the existence of an unstable localized eigenvalue (if this 
point lies in the region of consistent splitting), or else the presence of essential spectrum 
on the positive real axis. Therefore, given that D'{0) < a conclusion can always be 
made concerning the stability of the underlying solitary wave: a time-exponential instability 
to localized perturbations or else a convective-time-oscillatory instability (as described in 
the introduction and in Section ^ below) . It turns out that the sign of this derivative is 
intimately related to the underlying geometry of the ODE flow induced by the profile ODE. 
Indeed, for values of {c,q) such that there is a smooth saddle equilibrium at (ro(c, g),0) 
in the traveling wave ODE, let U~^{c,q;x) parameterize the stable manifold at (to,0) and 
U^{c,q;x) the unstable manifold. Define the Melnikov separation function 



(3.7) d{c,q) = det ( ,\, ^^+ ^_^, ) \x=o, 



Tx T~^ — T 

{fxY (r+-f- 

and notice that this represents a signed distance of U^ from C/+ along a normal section at 
U\x=o oriented in the direction of U^\x=o with respect to the right hand rule; see Figure pi 
As seen in the next lemma, the sign of the quantity -D'(O) is determined precisely by whether 
d(c, g) is an increasing or decreasing function of the wave speed c. 



Lemma 3.10. If d{c,q) = 0, then for the linearized St. Venant eigenvalue equations (3.1) 



about the corresponding homoclinic orbit of (2.4) we have sgnD'(O) = —sgndcd{c. 



Proof. This follows by a computation similar to those of |GZ| IZ21 IGoj . First, notice that 
since A = is a simple translational eigenvalue, we have -D(O) = and we may take without 
loss of generality W^ = W^ = U' at X = 0. Furthermore, a matrix perturbation argument 
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Figure 3: Depiction of the geometric interpretation of the Mehiikov separation function 
d{c,q), in the special case where Ux{^) = (0, 1)^- Here, note that, by a direct calculation, 
d{c,q) < 0. 



at A = shows the unstable subspace of the asymptotic matrix Aq(X) can be analytically 
extended to a neighborhood of the origin, with limiting values at A = given by the direct 
sum of the unstable subspace of Aq{0) and the vector (l, —g-riuo, TQ)gu{uo, to)~^,0) , with 
g as in (2.6) It follows then that we can choose the vector VF2~(";0) to be asymptotically 
constant with 



V^ = lim W^{x,0) = 1, 



-{s + l)uo 



,0 



(l,-(i/*(To),0)^ 



where the last equality follows from the relation Tq"*" Uq = 1. Moreover, assuming the 
homoclinic connection goes clockwise in the (r, r') variables and that //+ < and /u_ > 
denote the decaying modes of the asymptotic matrix j4o(0) at ±00, respectively, we have 



yf(o) 



-ro,cro,C;U_To 



-l^'j^ 



where we have used the fact that u = —ct for the vectors W^{-, 0) and W.^{-, 0). Therefore, 



det{V,-,V+,V, 



3 JIA=0 



(3.8) 



/ 


-To 1 To 


det 


cro -4f*(ro) -ctq 




^ cf^-To^ -ch+Tq 


c{ii- 


-^+)'^*l c -dMro)) 


c{fi^ 


-^+)((i/*(To)-c). 



which is non-zero by assumption. 

Continuing, since det(y]~, V^ , Vr^) is non- vanishing and det(W^f , W2 , W3'')|^_o,A=o = 
we have 

^ ^^' " A u\r- M- T/+^ 1^=0 



3 ANALYTICAL STABILITY FRAMEWORK 18 

and hence, to complete the proof then, we must calculate the numerator. Using the Leibnitz 
rule, we immediately find that 
(3.9) 
dx det{W^,Wi,W+)\x=o = det{dxWr,W^,W+)\x=o + ■■■ + det(Tyf , W^,dxW+)\x=o 

= det(C7., W^,dxW+ - dxWr)\x=o, 

where the vectors dxWj\x=o satisfy at 

(3.10) dxW'j\x=o = A{x,0)dxWj\x=o + U' 

for j = 1,3. Furthermore, differentiating the stable and unstable manifold solutions of the 
traveling-wave ODE at the endstate {tq,uo) (chosen in equilibrium) with respect to the 
wave speed c, we find 

/ d^f- \ 
dxWr\x=o = -\ dcu~ , dxW+\x=Q = - 

Note this is a reflection of the similar role of c and —A in the structure of the linearized 
traveling wave equation and in the linearized eigenvalue equation, respectively. 
Next, setting A = in (|3.1|) and integrating from x = — ootox = implies 




(n + cr) |a;=o = c - d/*(ro) 
for {t,u) = ((W2~)i' (^2~)2)) the first two components of the vector W2~(-,0). Similarly, 



recalling (3.10) it follows that 

{u + cr) |a;=o = / {u^ + cf^) dx = {u + cf) |°=_tc 

J±oo 



for (t, n) = {t ,u ). Therefore, returning to (3.9), we find that adding - times the second 

row to the first yields 

(3.11) 

/ i(c-d/.(ro)) ^0 \ 

dxdet{W^ ,W2 ,Wj^)\^=o,x=o= \ u^ * dc{u~ - u+) U=o 

\f-\u^y * f~^dc{u--u+yj 

= - (dMro) - c) det [--2^- J ^-2Q^^^- _ ^+y) U=o 

= -cf-2(0)((i/.(ro) - c) det (^ J^, ^^|_^+ ^ .^I j,) U.=o 
= -cf-2(0)((i/,(ro)-c)9,(i(c,5). 
The proof is now complete by dividing ( 3.11| ) by ( |3.8[ ) and noting that (^_ — ;U_|_) > 0. D 
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Corollary 3.11. sgndcd{c,q) < is necessary for stability of point spectrum of the corre- 
sponding homoclinic, corresponding to an even number of roots with positive real part. 



Remark 3.12. In the stable case considered in |N2| . for which periodics appear as c is 
decreased from the homocUnic speed Chom) we see that periodics exist in the situation that 
d < 0, which should be visible from the phase portrait. That is, the lower orbit coming 
from the saddle point should pass outside the upper branch. Likewise, we can check the 
sign for Jin-Xin by looking at the phase portrait for nearby periodic case. In fact, we can 



compute directly that dgd > for the Jin-Xin case, as done in Section [3^ below, yielding 
instability for that model. 

Remark 3.13. Note that dcd < implies for c < Chom! where Chom denotes the homoclinic 
wavespeed, that the unstable manifold at the saddle equilibrium spirals inward, either 
toward an attracting equilibrium in which case periodics don't exist, or toward an attracting 
periodic cycle, in which case the enclosed equilibrium is a repellor. Likewise, dcd > implies 
for c < Chom that the stable manifold at the saddle equilibrium spirals inward in backward 
X either toward a repelling equilibrium in which case periodics don't exist, or toward a 
repelling periodic cycle, in which case the enclosed equilibrium is an attractor in forward x. 
That is, generically, the stability condition dcd < is associated with the property that the 
enclosed equilibrium be a repellor, so long as decrease in speed from Chom is associated with 
existence of periodics. (If, rather, increase in speed is associated with existence of periodics, 
then we reach the reverse conclusion that Evans stability correspods to the property that 



the enclosed equilibrium be an attractor.) By Remark 2.4 therefore, we expect stability 



(of point spectra) of homoclinics in the St. Venant case studied in |N2] . 



As described in Remark 3.13, it is a curious fact that stability in temporal variable t of 



the homoclinic as a solution of PDE (2.2) is associated with stability in spatial variable x 
of nearby periodics as solutions of the traveling- wave ODE (2.4), and, therebyPJ instability 



of the equilibrium contained within the periodics. The latter is also the condition that 
the homoclinic itself be stable in forward j;-evolution as a solution of the traveling-wave 
ODE. It is not clear whether this is a coincidence, or plays a deeper role as a mechanism 
for stability. However, once we understand the existence theory (itself nontrivial, but often 
observed numerically), we may make conclusions readily about stability just by examining 
stability or instability of the enclosed equilibrium as a solution of the traveling- wave ODE. 

3.4 The viscous Jin-Xin equations 

While of theoretical interest, unfortunately the stability index dcd derived above appears 
to be analytically incomputable for the general St. Venant equations considered here (al- 
though, one can use the "rule of thumb" described above). Next, we consider as simpler 



'^Assuming the generic case where there is only one limit cycle generated from increasing or decreasing 
the wave speed c, with an associated return map that is strictly attracting or strictly repelling. 
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related example the viscous Jin-Xin model 



(3.12) 



Tt -U^ 



0, 



ut - c^Tx = -fir) -U + Uxx 



for which we can directly calculate the value of dcd: as we will see, Lemma 3.11 holds in 
this simpler setting as well. Considering traveling waves in the moving coordinate frame 
X — ct gives traveling-wave ODE 

{c^-ci)T' = -fiT)+cr-q-CT", 

where u = q — ct. Taking c = Cs, this equation is seen to be Hamiltonian and reduces to 
the nonlinear oscillator 



(3.13) 



CT 



-/(r) + CT-q, 



for which all periodics and the bounding homoclinic arise at once. Indeed, this equation 
may be integrated via quadrature: multiplying by r' and integrating, we obtain 



(3.14) 



(c/2)(r 



/\2 



-F(r) + CT^/2 -qT + H, 



F := J f, which is of standard Hamiltonian form with ODE energy H and effective potential 
energy V{t; c, q) = F{t) — |r^ + gr. In particular, we see that the homoclinic wavespeed 
Chom is precisely Cg given in (3.12). 



Denoting the roots of the equation V'{t) = as tj, we have that orbits through the 
equilibrium solutions (tj,0) of (3.14) line on the level sets V{t) = H, and correspond 



to homoclinic orbits provided that V"{tj) < 0. In this case, the orbits of (3.14) can be 
implicitly defined via the relation 



(3.15) 



r' = ±V(2/c)(-F(r) + ct^I -qT + H) 



from which solutions of (3.12) may be obtained solving (3.15) for r as a function of x 
and using that u = q — ct. Furthermore, a direct calculation shows that the sign of 
V"{tj) = —df{Tj) + c alternates as Tj increases, and hence we find that the subcharacteristic 
condition is satisfied if and only if V"{tj) < 0. It follows that all non-trivial solitary wave 
solutions of (3.12) are stable to low-frequency perturbations, and hence have stable essential 
spectrum. 

The analysis of the point spectrum of the corresponding homoclinic orbits is handled by 



the usual Evans function framework as described in Section 3.2 To this end, notice that in 
co-moving coordinates the linearized equations are 



(3.16) 



n 



CTx 



0, 



Ut - CUx - cItx = -df{T)T -U + Uxx, 
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where U = (r, u) denotes background profile. Taking the Laplace transform in time imme- 
diately leads to the linearized eigenvalue ODE 



(3.17) 



AT — CT — U 
2 / 



0, 



\ / 2 / 

Xu — CU " ~ 



-df{T)T — U + u" , 



governing the spectral stability of the underlying wave to localized perturbations. In order 
to use the abstract Evans function setup of Section |3.2[ we write the eigenvalue system 



(3.17) as a first-order system of the form W' = A{x, \)W , where 
(3.18) 



W :-- 



A:-- 



\/c -1/c 

1 

cl\/c + df{f) A + 1 cl/c-c/ 



In order to apply Lemma 3.4, and hence define the corresponding Evans function D{\) as 



in (3.5), we must verify that we have extended consistent splitting in this case. 



To obtain consistent splitting in the homoclinic case, we must have stability of constant 
solutions at the endpoint r = from which the homoclinic originates, which as noted earlier 
is equivalent to the subcharacteristic condition 

df*{Tj) < Cs, 

where /*(t) := —u^{t) = /(r) is determined by the value u*(t) = —f{u) at equilibrium, 
hence df^.{T) = df^r). Thus, the subcharacteristic condition is df{T) < Cg, and so as 
discussed above it is satisfied at alternate roots. (Recall that, by our previous analysis, 
instability of the constant solution is necessary for Hopf bifurcation and thus for existence 
of a nonlinear center and bounding homoclinic.) With this setup, we see that Lemma 



3.4 applies to the eigenvalue system corresponding to (3.17) and hence we may define the 



Evans function D{X) as in (3.5). Furthermore, a direct calculation shows that the proof 
of Corollary 3.11| carries over line by line to the viscous Jin-Xin case considered here. It 
follows then that we can determine the stability of the solitary wave by direct considerations 
of the Melnikov separation function d{c, i 
geometric nature of the function d{c, 
at the homoclinic. 



defined in (3.7). It turns out that due to the 



the sign of the derivative dcd is trivial to compute 
Indeed, note that for general c the traveling-wave ODE is 



(3.19) 



CT 



-/(r) + CT- q- {c^ 



2\ / 

cJt . 



It follows that H' = — (c — c^)|r'| , where H is the the Hamiltonian defined as in (3.14), and 
hence, for c > Cg, the ODE flow decreases the Hamiltonian H. Noting that H is minimized 
at the equilibrium enclosed by the Homoclinic, it follows that we must have d > for 
< c — Cs <^ 1 and hence we conclude that dcd > 0. By Corollary 3.11 then, we conclude 



that the solitary waves of (3.12) always have unstable point spectrum. 



This example, although much simpler than the full St. Venant equations, illustrates the 
power of being able to relate the stability of the solitary wave as a solution of the governing 



4 HIGH-FREQUENCY ASYMPTOTICS 22 

PDE to information about Melnikov separation function d(c,q), an inherently geometric 
quantity encoding information about the ODE phase space. For the St. Venant equations, 
however, we must resort to numerics to study the point spectrum of the hnearized operator. 
This analysis is carred out in the following two sections. As we will see, however, although 
we are not able to analytically compute the stability index sgndcd{c, q) in this general case, 
our numerical experiments are consistent with the rule of thumb relating the PDE stability 
to the stability of the enclosed equilibrium solution of the ODE. 

4 High-frequency asymptotics 

Our next goal is to carry out a numerical Evans study of the point spectrum for the lin- 
earized equations ( |3.1| ) for homoclinic orbits corresponding to various choices of the turbu- 
lent friction parameters (r, s). Indeed, later we will see that a winding number argument 
allows us to numerically compute the number of roots of the Evans function contained in 
a given compact region of the complex plane. The goal of this section is thus to eliminate 
the possibility of arbitrarily large unstable eigenvalues, so as to reduce our computations 
to a compact domain in A |Brl IBrZl IHLZ| HLyZl HLyZ2| . This can be done by a num- 



ber of techniques, including energy estimates and high-frequency asymptotics. As energy 
estimates degrade with small viscosity, we opt to use the latter techniques, which is well 
adapted to the situation of multiple scales such as are involved in our examples jPZl IZ5j . 

4.1 Approximate block-diagonalization 



To begin, we show that the coefficient matrix ^(a;. A) in the Evans system (3.2) can be 
written, for |A| S> 1, as a block diagonal matrix plus an asymptotically small term. We 
begin by expanding W'{-) = A{-, X)W{-) a.s A = B + R, where 



B :-- 
(4.1) 

R:-- 



A 


0\ 



















f2 


: 














-aX 
cv 


\ o/ 
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/o 





k 

















= : 








(s+l)f={q-cf)'-- 


-OLx 


rr 


+Hq-cfr-' 


— cf^+af^/c 


/ 


m. 


n 



where a and ax are defined in (3.3) and where r" is then computed using equation (2.4) 



Since R = 0(1) as |A| — )• oo, we expect that the behavior of the Evans system (3.2) is 
governed by the principle part B(x, A) for A sufficiently large. Notice, however, that it is 
not straightforward to characterize this due to the multiple spectral scales associated with 
the growth of the eigenvalues of B. Indeed, a straightforward computation shows that the 

eigenvalues of the principle matrix B are given by - and if a/-, and hence the spectrum 



of B has two principle growth rates: order A and order A^'^. In order to keep track of both 
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of these scales, one must make a number of careful transformations, which we now describe 
in detail. 

First, we introduce the transformations 

'1 0\ / 1 0^ 

(4.2) T := ( 1 , T^i = 10 

1/ \-e Oy 

where 9 := — -. Then we readily see that 

/^00\ / ke Ok 

(4.3) T-^BT = lof"^ r^ , T-^RT = 

\ ^ 0/ \l + ne-ke'^ m n-kOj 

and, defining 

0\ /A/c 

(4.4) Bi:=T-^BT- [f^e 0=0 r^ | , W:=TWi, 

0/ V ^ 

we have W[ = {Bi + Ri)Wi, where 

/ 
Ri:=T-^RT-T-^T' + If^e 

\ 

(4.5) 

ke k \ /ki h 

f'^6 =: mi 

I + nO — kO"^ — 6' m n — kOj \ni oi pi^ 

Now, define 

/I \ /I 

(4.6) Ti := 1 , T^-i = 1 

\0 ;U^/A/ yo (Ai\/A)-i^ 



where u = ^^ and notice that 



tWi/ 



\ ki hfiVX 



(4.7) T^^BiTi =0 \/A^f2 , Tf ^i?iTi = | mi 

VO \/A^f2 / 

Setting Wi = T1W2, we thus have W^ = (^2 + ^2)W^2, where 



VA^f^ / V^ ^ ^1 



'A 



hfiVX' 



(4.8) ^2 := I y/Xfif 

,0 \/A;ur2 
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and 

/O hfiVX'^ 

(4.9) i?2 := T^^RiTi - T{^T[ - 

\0 



with p2 := Pi - ^ = Pi + V- 



ki 



mi 

rii oi 



0' 

01 



P2j 



Continuing, define 
(4.10) Ts := 




t; 




and set B^ := T2 .622^2 and W2 := T2W3. Then by a direct calculation we have 



(4.11) 
where 

(4.12) 

and 

(4.13) 



W!, = iB3 + R3)W3, 



B^ 



lijiy/X hfiyX 
/A^f2 

-VXfif^ 



/a 

c 

r 



A 
c V '^ 
A 



a\ 



P rp— 1 T-> rp rp—lrpi 

K3 ■— ±2 ^2J^2 - J^2 ^2 



k 







V 















mi 



+ni/A/A/j P2+oi/\/A^ P2-oi/'/Am 



2 2 2 

-mi+ni/>/A^ P2+oi/>/A/J P2^oi/>/Am 

2 2 2 

Finally, we introduce the transformations 
/l 



'A;3 
^3 ms "-3 I . 
,03 7713 71,3, 



(4.14) 



n-.-- 



1 
\o 



A VA 



1 



^-1 

^3 



Va 
1 



.0 



where 9 = ^ and note that 



T^^Bsn 



r ^^ 
\0 



-^/^/ 



(4.15) T^^R3n = 

(k3-{h + os)j^ (k3-2ms-{h + os)j^)j^ {ks-2ns-{k + o^)j^)j^\ 



V 



^3 "^3 + ^3;^ "3 + ^3^ 



03 



777.3 + 03 



Va 



m + 03 



yx 



/ 
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Defining tlien 



(4.16) 



and 



Ba 



(- 

I c 

r 



\ 







(4.17) 



Ra 



T^^RsT. + 








it follows that the original Evans system (3.2) can be written in the approximately block- 
diagonal system 



(4.18) 



W' = iB4ix,X) + R4ix,X))W 



where R4 = o(A^) as A — )• 00 for all e > 0. With this result in hand, we can now obtain 
bounds on the size of the unstable eigenvalues A corresponding to the linearized St. Venant 



equations (3.1). Before stating these bounds, however, we need the following abstract 



tracking lemma concerning the eigenspaces of approximately block-diagonal systems. 
4.2 Quantitative tracking lemma 



We now recall the following quantitative estimate |HLyZl[ |HLyZ2[ . Consider an asymptot- 
ically constant approximately block-diagonal system W' = (M -|- @){x)W, 



(4.19) 



W 



W4 



M 



M_ 

M4 



G 






with KM_ < c_/ and KAf+ > c+I satisfying 

(4.20) c+-C-> 5{x) > 0, 

where 3f?M := (l/2)(Af + M*) denotes the symmetric part of a matrix M. 
Denote by 



(4.21) 



C±{x) :-- 



10— \ 



le. 



2|e. 



ie_. 



le. 



219. 



19- 



19. 



the roots of 

(4.22) PiS,x) := ( - 5 + |9__| + |9++|)C + I9-+I + |9+_|C^ = 0, 

where, here and below, | • | denotes the i"^ matrix operator norm. 
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Lemma 4.1 ( |HLyZl[|HLyZ2| ). Suppose that 

(4.23) (^> |e__| + |G++|+2^eI 



le. 



Then, (i) < C- < C+; (''-'i') ihe invariant subspaces of the limiting coefficient matrices 
(M + 0)(iboo) are contained in distinct cones 



n- 



{\W-\/\w+\ < c-}, n+ = {|T^_|/|t^+| > c+}, 



and (Hi) denoting by S^ the total eigenspace of (M + 0)(+oo) contained in r2+(+cxD) and 
U^ the total eigenspace o/(M + 0)(— oo) contained in $7_(— oo), the manifolds of solutions 
of (4.19) asymptotic to S~^ at x = +oo U~ at x = — cxd are separated for all x, lying in 0+ 



and 0_ respectively. In particular, there exist no solutions of (4.19) asymptotic to 5"*" at 
+00 and to U~ 



at 



-oo. 



Remark 4.2. In the case that c_ < at +00 and c+ > at x = — 00, S~^ and U^ 
correspond to the stable subspace at +00 and the unstable subspace at —00 of the limiting 
constant-coefficient matrices, and we may conclude nonexistence of decaying solutions of 



(4.19), or nonvanishing of the associated Evans function. More generally, even if (as here), 



the property of decay is lost as we traverse the boundary of consistent splitting, so long as 



there remains a positive spectral gap 6, and (4.23) remains satisfied, we may still conclude 



nonvanishing of the Evans function defined as above by continuous extension of subspaces 
U- and S+. 



Proof ( |HLyZli[HLyZ2l ) From ^J^, we obtain readily 
(4.24) 



\w^\' <c_|Ty_| + |e, 
\w+\' > c4-\w+\ - \e. 



\W- 



+ |e_ 






from which, defining (^ := \W- 
Riccati equation (' < P{(^,x). 



\W^\, we obtain by a straightforward computation the 



Consulting (4.22), we see that C < on the interval ("_ < (^ < ^+, whence O. 



{C < C_} is an invariant region under the forward flow of (4.19); moreover, this region is 
exponentially attracting for ( < (^+. A symmetric argument yields that r2+ := {^ > ^+} 



is invariant under the backward flow of (4.19), and exponentially attracting for C > ^_. 
Specializing these observations to the constant-coefficient limiting systems at x = —00 and 
X = -|-c«, we find that the invariant subspaces of the limiting coefficient matrices must lie 
in r2_ or r2_|_. (This is immediate in the diagonalizable case; the general case follows by a 
limiting argument.) 

By forward (resp. backward) invariance of $7_ (resp. i^+), under the full, variable- 
coefficient flow, we thus find that the manifold of solutions initiated along U^ at x = —00 
lies for all x in r2_|_ while the manifold of solutions initiated in S^ at x = -|-oo lies for all x 
in Q-. Since 0_ and 0+ are distinct, we may conclude that under condition (4.23) there 
are no solutions asymptotic to both U^ and 5"+. D 
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4.3 Bounds on unstable eigenvalues 

With the above tracking lemina in hand, we are able to state our main result concerning the 



confinement of the unstable eigenvalues of (3.1 ). After a change of coordinates switching 



and +, using notations of Section 4.1, equation (4.18) determines a system of form (4.19) 
with 



(4.25) 



G+4 



f 




M_ 



+ A^^/^ei , + A-^el , + A-3/20 



mi g2 

V 2 2 



+ 



A-1/2 




-e' + 9{ki-p2) 



v-l 







+ A-3/2 



2/^ 






2 
22 ni 





e+_ = G° + A-i/201 ^ A-^el_ + A-3/203 



+ 



A-i/2 



+ 

2/^ ^ 



P2j 



i rrti 
2 







e_- 



e. 



~ \ 2 2 / "I" '^ \^ 2/i 2^4 

e^_ + A"^/2e^_ + A"^e^_ 
= (f)+A-^/2f-^-^^')+A-i 



mi 
2 



£1. 
2/i 



+ A 



ini 
2^ 



-1 



+ A-3/2 



a2ni 





2/x 



for which evidently 6 > tK-v/- > |A|^/^— £= when 3f?A > 0. Then a direct application of 



Lemma |4. 1| yields the following high-frequency bound. 



Corollary 4.3. There are no unstable roots A, with KA > 0, of D with |A| > max^^ R^, where 
R is the only positive root of X^ — oqX^ — ai/2X^ — aiX^ — 03/2^^ — 02^^ — ^5/4^ ~ 0,3 
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with 



(4.26) 



ao 



H/2 



ai :-- 



J3/2 



a2 ■■-- 



O5/4 : = 



03 : = 



\/2i;(|eo_ 


+ 




+ 2y 


^leo+iie^ 


\ 
/ 
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^2^2^90 


-+ 


©V-l 


+ |ei 


+l|0°+-l 





2^ |Gi_| + |ev+| + 2 J|eL+||e^_| + |eo_+||e^_| + le^ ^||eo.. 



^2^2^90 


+II0M 


+ |e^ 


le^ 


_l + 


e^ 


-liev 


1 






1 + 10^+ 


1 + 2^ 












V2^(|G2_ 


eM 


le^ 


_| + 


eu\ 


let 





2z.2,/|e2 1163 I 



2z/|Gi]_+| 



In particular there are no unstable roots X with modulus larger than W' where R is the 

only positive root of X^ — (maxao)X^ — (niaxax/2)^^ — (maxai)X'^ — {iaay.a^i2)X^ — 
^2 



(maxa2)X — (max 05 m)^ — max 03. 



To aid in the utilization of the bounds in Corollary 4.3, note that R is smaller than 



• the square of the only positive root of 
(4.27) 



=: X^ - aoX^ - aiX^ - agX - 03; 

• the only positive root of 

(4.28) X"^ -{al + 2ai + a2)X-{a2 + 2a^). 

Both bounds can be explicitly given. Here, however, we determine R by numerical solution 



of the full eighth-order polynomial equation given in Corollary (4.3). 



Proof. By Lemma 4.1 , the triangle inequality and subadditivity of the square-root function, 



5 NUMERICAL INVESTIGATIONS 



29 



there are no unstable eigenvalues A with 
(4.29) 



|A|i/2 > 



^(|e^_| + |e%|+2,/|e^^liel 



+ 


A| 


+ 


A| 


+ 


A| 


+ 


A| 


+ 


A| 


+ 


A| 



-l/4_ 



V2^2J\@o \\@i \ + |e^ 1190 I 



_^^^V2^{\el_\ + \QU\ + 2^|eL+||ei._| + |G0_+||e2._| + |e^+||e^. 



-3/4_ 



\/2^2 Jieo 1163 I + leL , 1162 I + 162 , 1161 I 



^ V2^ {\el_\ + |6^+| + 2^|e2_+||6i_| + |6i_,||6i 



-5/4 



-3/2; 



2i.2J|62,||63 I 



/2l^ie3 



also written |A|i/2 > a^ + ay2\M~^^^ + ai\X\-'^^^ + a3/2\M~^^^ + a2|A|-i + a5/4|A|-5/4 + 
ffl3|A|~3'2. rpj^g condition is satisfied as soon as |Ap^ > max^: i?. D 

For those interested in reproducing or using the high frequency bounds obtained above, 
we record in the Appendix the formulas for the various terms listed above in terms of the 
original underlying homoclinic profile. 



Remark 4.4. By similar estimates (see the abstract Tracking Lemm,a of JMaZ3[ \PZ[ \Z3^ ). 
one may obtain as in \HLyZl^ the asymptotic description 



(4.30) 



D{\) ~ cie'^^VA 



for some ci,C2 G M as |A| — t- 00, with convergence at rate 0{\\\ ^''^). 



5 Numerical investigations 

In this section, we carry out a numerical study of the stability of the solitary wave solutions 
of the generalized St. Venant equations. In particular, we perform numerical Evans func- 
tion calculations to detect the presence of unstable point spectrum. In the two examples 
considered, we find that one has stable essential spectrum with unstable point spectrum, 
while the other has stable point spectrum with unstable essential spectrum. Furthermore, 
for the latter case we perform a time evolution study to illustrate the convective nature of 
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the instability arising from the essential spectrum. For illustrative purposes, however, we 
choose to perform a numerical Evans study of the solitary wave solutions of the viscous Jin- 



Xin equations. Recall that by our considerations in Section 3.4, we expect the linearization 
to have unstable point spectrum, in particular, a positive real eigenvalue, corresponding to 
exponential instability of the solitary wave. 

5.1 Jin— Xin Example 

To begin, consider first the set of equations 

Tt - Ux = 0, 



Ut- Tx = T ^/^ - n + Uj 



(5.1) 

corresponding to (3.12) with nonlinearity /(r) = — r^^'^, mimicking the St. Venant struc- 



ture, and Cs = 1. For this example, we consider the solitary wave profile U = {t,u) satisfying 
the ODE 

l//\2 Iq lo 1 

2 ^ '' 3 2 2 



corresponding to the profile ODE (|3.14) with {q,H) = (2,-2). This profile is numerically 



generated using MATLAB's boundary- value solver bvp5c |KLj and is depicted in phase 
space Figure [4]j^a) and is plotted against the spatial variable x in Figure [4][^b). To study the 
spectral stability of this profile, we must next consider the linearized eigenvalue equations 



(3.17) with Cs = 1 and df{T) = — 2r. Since we expect the existence of a real unstable 



eigenvalue by the calculations in Section 3.4, we numerically compute the Evans function 



along the positive real line and show the existence of a positive root. This is accomplished 
here by using the polar-coordinate method of |HuZj as well as other standard procedures; 
see also [BrZj . In particular, the Evans function was computed in MATLAB using RK45 
to solve the associated ODE's. The adaptive error control provided by RK45 gives an 
estimate on truncation error which, by the results of |Z6] . then translate to convergence 



error estimates. Furthermore, we find analytically varying initializing bases {R- } for the 
Evans function computation via the method of Kato; see \G7j\ IHuZl IBrZllBHZ] . Following 
this procedure then, we plot in Figure |4] (c) the output as evaluated on the interval [0, 1] 
and clearly see that D{X) vanishes (geometrically) twice within this interval: once at A = 
corresponding to the simple translational eigenvalue and once at a positive Aq > yielding 
the expected instability. 

An alternate way of numerically seeing this instability through the use of an Evans 
function is through the use of a winding number calculation. Indeed, the analyticity of the 
Evans function on the spectral parameter A implies the number of solutions of the equation 
D[\) = within the bounded component of a contour F, along which D|r is non-vanishing, 
can be calculated via the winding number 

^ ' 2m Jr -D(A) 
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{ay 




(b) 




(c) 




Figure 4: Evans function output for the Jin-Xin equations with /(r) = — r^^'^ and Cg = 
1, q = 2, and H = —1/2. We display the honiochnic profile in phase space in Figure 
(b) and against the spatial variable in Figure (b). In Figure (c), we observe the Evans 
function output evaluated on the real line, illustrating the existence of an unstable positive 
eigenvalue. 



Notice that unlike our previous Evans function calculation, which detected only real roots of 
D{X) = 0, this method allows one to exclude the possibility of unstable eigenvalues within 
a given compact domain. This advantage will be exploited below in conjunction with the 



eigenvalue bounds of Corollary 4.3 to conclude not only the absence of real unstable eigen- 



values but of any unstable eigenvalues. To illustrate the numerical procedure, we consider 
again the Jin-Xin profile studied above. Notice that finding the analytically varying initial- 
izing bases {i? • } as described above preserves analyticity of the Evans function allowing us 
to use, via the argument principle, a winding number computation to determine the pres- 
ence of roots inside a given contour. For this example, we choose a semi-circular contour of 
radius one and inner circle radius 10~^. Notice for a winding number calculation the contour 
must not go through the origin since D{0) = 0, and hence we must deform the standard 
semi-circular contour slightly near the origin. Using sufficient number of mesh points then, 
we can ensure the relative change in D{Xj) between consecutive mesh points {Aj} is less 
than 0.2, thus assuring an accurate winding number count by Rouche's Theorem. For this 
specific example, a plot of the Evans function output restricted to the semi-circular contour 
described above is shown in Figure [5J From this, it is easily seen that the winding number 
is one, again verifying the expected instability. 

Next, we perform the analogous numerical study for the St. Venant equations (2.2) 



for two distinct values of the turbulent friction parameters. First, we consider the case 
(r, s) = (1,1) demonstrating stable essential spectrum and unstable point spectrum. Second, 
we consider the common case (r, s) = (2, 0) demonstrate unstable essential spectrum and 
stable point spectrum. As mentioned in the introduction, this instability is convective in 
nature and we demonstrate this with a time evolution study. 
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(a) 





ib) 



(c) 



Figure 5: Evans function output for the Jin-Xin equations with /(r) = — r~^" and Cg = 1, 
q = 2, and H = —1/2. Figure (a) depicts our semi-circular contour used for the winding 
number calculation. In Figure (b), we have the Evans function output evaluated along the 
contour and in Figure (c) we zoom in on the origin. The winding number is one for the 
contour displayed, as expected. 



5.2 Case (r,s) = (1,1) 
First, consider the equation 



(5.3) 



ut + ((0.444)-V 



2\ 



0, 

1 - T^ii + 20(r" 



corresponding to the St. Venant system (2.2) with (r, s) = (1, 1), F = 0.222, and v = 20. 
In this example, we consider the traveling wave profile associated with the corresponding 
traveling wave ODE (2.4) with wave speed c ~ 3.1869 and integration constant q = 1 + c. 



Following the protocol previously described, in Figure ^^a) we depict the corresponding 
orbit in phase space along with the Evans function output along the real line. From Figure 
p[b), it is clear there exists an unstable real eigenvalue associated with this profile. Al- 
though not depicted here, this was further verified by a winding number calculation as in 
the Jin-Xin case above, in which one finds a winding number of one for an appropriate 
semicircular contour. Moreover, notice that the limiting constant state (1,0) satisfies the 
subcharacteristic condition 

3/2 

2VF < ( lim t(x) ' 

corresponding to spectral stability of the limiting state. It follows that although this profile 
admits unstable point spectrum, its essential spectrum is stable. We verify this numerically 
using the SpectrUW package developed at the University of Washington [DK| . which is 
designed to find the essential spectrum of linear operators with periodic coefficients by using 
Fourier-Bloch decompositions and Galerkin truncation; see [CuD|, ICDKKl IDKj for further 
information and for details concerning convergence. Although this package is designed 
for periodic problems, the essential spectrum corresponding to the limiting constant state 
(treated as a periodic function) of our homoclinic is precisely the essential spectrum of the 
original homoclinic profile, and hence we can numerically compute the essential spectrum 
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of a solitary wave using SpectrUW in this way. In particular, we find that the essential 
spectrum is stable, as expected; see Figure [gFc). 



(a) 




(b) 



(c) 



Figure 6: Demonstration of unstable point spectrum in generalized St. Venant system. 
Here r = s = 1, u = 20, F = 0.222, c ~ 3.1869, and q = 1 + c. We have the phase portrait 
plotting T verse r' in Figure (a), and the Evans function evaluated along the real line in 
Figures (b). Finally, Figure (c) depicts the numerical evaluation of the essential spectrum 
of the homoclinic using the SpectrUW package. 

We now wish to connect this numerical experiment to the analytical Evans function 
calculations of Section [3. 3[ in particular Corollary |3. 11 and Remark 3.13 To this end, we 
note that numerically it is found that periodics exist for the system (5.3) as c is increased 
through the homoclinic wave speed, which here we have chosen as Chom ~ 3.1869. Further- 
more, as seen by Figure M^a) the equilibrium solution enclosed by the homoclinic orbit is 
a repeller. Together these considerations suggest that dcd{c, q) > hy our discussion in 
Remark 3.13 Note this is verified by noting from Figure p[b) that -D'(O) < and using 
Lemma 3.10 Thus, the instability of this profile follows trivially from direct analysis of the 
corresponding phase portrait and Corollary |3.11[ Notice, however, that this example shows 
the dependence on the "repeller/attractor" rule of thumb on how the periodic orbits are 
generated from the homoclinic. 



5.3 Case (r,s) = (2,0) 
Next, we consider the equation 



(5.4) 



Wi + ((18)-V 



2\ 



0, 
1 



ru^ + 0.1 (r^^ii. 



corresponding to the St. Venant system with (r, s) = (2,0), -F = 9, and z^ = 0.1. Here, we 
consider the solitary wave profile associated with the corresponding profile ODE with wave 
speed c ~ 0.7849 and integration constant g = 1 + c, which is depicted in phase space in 
Figure ma). We begin by discussing the point spectrum associated to this solitary wave 
solution. Noting that periodics are seen to be numerically generated as the wave speed is 
decreased from the homoclinic speed, the fact that the equilibrium enclosed in the homo- 
clinic is a repeller suggests that dcd{c,q) < by Remark 3.13 In particular, our heuristic 
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suggests that the stability index derived in Section 3.3 does not yield any information con- 
cerning the existence of unstable point spectrum of the corresponding linearized operator. 
This motivation is further supported noting from the Evans function output along the real 
line in Figure [7[b) that -D'(O) > 0, and hence sgndcd{c,q) < by Lemma 3.10 as expected. 
In order to further study the point spectrum of the corresponding linearized operator, we 



need a high frequency bound on the eigenvalues of (|3.1|). By using the numerically computed 

implies that any 



4.3 



profile to evaluate the functions aj in (4.26), an application of Corollary 
unstable roots of D(X), if they exist for the given parameter values, must satisfy the bound 
|A| < 308. In FigureJTJ^c) and (d) the output of the Evans function evaluated on a semi-circle 
of outer radius 308 and inner radius 10~^, from which we compute that the winding number 
(5.2) is equal to zercrl In particular, it follows that the linearized eigenvalue problem 



in consideration has stable point spectrum as expected. As noted in the introduction, 
this example is notable as being the first example of a second-order hyperbolic-parabolic 
conservation or balance law which admits a solitary wave with stable point spectrum. 



Remark 5.1. The eigenvalue estimates of Corollary 4-3 should be contrasted with the results 



of a convergence study based on (4.30). In the example considered here, the results of a high 



frequency convergence study are listed below, with "Relative Error" denoting the maximum 



relative error between the numerically computed value of D{X) and an approximant (4.30) 
determined by curve fitting on the semicircle |A| = R, KA > 0.' 



itive Error 


Radius R 


0.1112 


8192 


0.1525 


4096 


0.2062 


2048 


0.2740 


1020 


0.4230 


308 


0.4522 


250 


0.5484 


130 



The radius R = 308 derived from our high frequency asymptotics and tracking in Corollary 
thus correspond to a relative error of approximately R = 0.4230, and hence the estimates 



in Corollary 4-3 are quite efficient even though they produce a radius that is relatively 
large. Furthermore, from this convergence study it seems then that the convergence rate of 
OdAI"^'^) from (4.30) breaks down for R ss 250. This suggests that the bound of R = 308 
obtained from our tracking estimates in Corollary 4-3 is very close to the boundary of the 
high frequency asymptotics. 

Next, recalling that the subcharacteristic condition reduces to F < 4 in this case we see 
that the solitary wave must have unstable essential spectrum. This is verified numerically 
in Figure me), where again we approximated the solitary wave with a periodically extended 



®In this example, the Evans function was evaluated on the contour with 3457 mesh points, yielding a 
maximum relative error between mesh points of 0.2. As stated before, this is sufficient to yield an accurate 
winding number by Rouche's Theorem. 
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version of itself with very large period and used package SpectrUW. However, by perform- 
ing a time evolution study of the generated homoclinic profile we find that the nature of 
this instability is seen to be convective. That is, a small perturbation leaves the original 
profile relatively unchanged in shape, but grows as an oscillatory time-exponentially grow- 
ing Gaussian wave packet after emerging from the profile; see Figure [8J For this numerical 
study, we used a Crank-Nicholson finite difference scheme, i.e. we used a forward difference 
time derivative and a centered, averaged spatial derivative approximation. This solitary 
wave is thus seen to be metastable in the sense that the perturbed solitary wave propagates 
relatively unchanged down the ramp, while shedding oscillatory instabilities in its wake. In 
the next section, we further examine this phenomenon and its relation to the stability of 
nearby periodic waves. 



(a) 





ib) 




id) 




(e) 




Figure 7: Evans function ouput for St. Venant system with r = 2 and s = corresponding 
to a homoclinic profile. Here c ~ 0.7849, q = 1 + c, F = 9, and i/ = 0.1. In Figure (a) 
we have the phase portrait plotting r verse r'. In Figure (b) we plot the Evans function 
evaluated on the real line. In Figures (c) we show the Evans function output evaluated on 
a semicircle of radius R = 308 with a small inner circle of radius 10~^, where we zoom in 
on the origin in Figure (d). The winding number is zero in these computations. Finally, 
Figure (e) depicts our approximation of the essential spectrum as generated by SepctrUW. 



6 Dynamical stability and stabilization of unstable waves 

As noted in the introduction, similar metastable phenomena to that observed in Section 



5.3 have been observed by Pego, Schneider, and Uecker |PSUj for the related fourth-order 
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(a) 



(b) 



(c) 



(d) 




Figure 8: Time evolution snap shots for St. Venant equations with r = 2, s = 0. In 
Figure (a) we have a perturbation to the right of the profile solution. In figure (b) we see 
the perturbation moving left emerge from the profile and grow as Gaussian curves displayed 
in Figures (c) and (d). The solution settles to a translate of the profiles (thin lines). Here 
c ^ 0.7849, fl = 1 + c, F = 9, and z^ = 0.1. 



diffusive Kuramoto-Sivashinsky model 



(6.1) 



ut + d^u + d,j.u + 



drU^ 



0, 



an alternative model for thin film flow down a ramp. They describe asymptotic behavior of 
solutions of this model as dominated by trains of solitary pulses. Indeed, such wave trains 
are easily observed on any rainy day in runoff down a rough asphalt gutter as found by 
many roads in the U.S. (e.g., near the last author's home), as are small oscillations between 
pulses as might be suggested by the analyses here and in |PSUj , corresponding to convective 
instabilities. 

Thus, both of these sets of results may be regarded as partial explanation of the some- 
what surprising phenomenon that asymptotic behavior of actual inclined thin-film fiow 
is dominated by trains of pulse solutions which are themselves unstable. Regarding this 
larger issue, we have a further observation that we believe completes the explanation of this 
interesting puzzle, at least at a level of heuristic understanding. 

The key is to reformulate the question as: how can we explain observed stable behavior 
of trains of solitary waves that are in isolation exponentially unstable? Or, more pointedly: 
how can a train of solitary pulses stabilize the convective instabilities shed from their neigh- 
bors? Phrased in this way, the question essentially answers itself: it must be that the local 
dynamics of the waves are such that convected perturbations are diminished as they cross 
each solitary pulse, counterbalancing the growth experienced as they traverse the interval 
between pulses, on which they behave as perturbations of an unstable constant solution. 
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This diminishing effect is clearly apparent visually in time-evolution studies; see Figure |8| 
However, it is a challenge to quantify it mathematically. In particular, it is only partially 
but not wholly encoded by the point spectrum that we usually think of as determining local 
dynamics of the wave |Hel IGZ| IZHj . 

Rather, the relevant entity appears to be the dynamic spectrum, defined as the spectrum 
of the periodic-coefficient linearized operator about the periodic wave obtained by pasting 
together copies of a suitably truncated solitary pulse. Here, the choice of truncation is not 
uniquely specified, but should intuitively be at a point where the wave profile has "almost 
converged" to its limiting endstate. This dynamic spectrum would govern the behavior of 
an arbitrarily closely spaced array of solitary pulses, so captures the diminishing property if 
there is one. A bit of thought reveals the difficulty of trying to capture diminishing instead 
by decrease in some specified norm. For, the decay we are expecting is the diffusive decay 
of a solution of a heat equation, which does not occur at any specified rate when considered 
from a given norm to itself, but shows up in long-time averaged behavior from a more 
localized (e.g., L^) norm to a less localized (e.g., L^ or L°°) norm, and whose progress is 
difficult to measure by a "snapshot" at an intermediate stage. 

In Figure [9] we depict the dynamic spectrum of the profile used in the time evolution 
study of Section |5.3[ along with the periodically extended version of this profile used to 
carry out the necessary numerics. We see that the dynamic spectrum indeed appears to 
be quite stable, despite the instability of the essential spectrum of the linearized operator 
about the wave (corresponding to the convective instabilities discussed above), supporting 
our picture of pulse profiles as "de- amplifiers" that can mutually stabilize each other when 
placed in a sufficiently closely spaced wave train. 



(«) 



ib) 
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Figure 9: In Figure (a), we plot the periodic spectrum of the periodically extended version 
of the profile found in Section |5.3| given in Figure (b). This dynamic spectrum of the 
associated homoclinic profile is seen to be stable, in contrast to its essential spectrum which 
was seen to be unstable in Figure [7| This spectrum was found using again the SpectrUW 
package developed at University of Washington. Note that the dynamic spectrum includes 
also small (stable) loops close to the point spectrum of the homoclinic [G] . 



This suggests also that there should exist nearby periodic wave trains that are stable, 
despite the instability of a single solitary wave. In companion papers [B JNRZlj IB JNRZ2J . 
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we show that this is indeed the case, estabhshing spectral, hnearized, and nonUnear stabil- 
ity of periodic wave trains in a band near the homoclinic limit. We point to the dynamic 
spectrum, or similar quantification of de-amplification properties, as an interesting direction 
for further development; see [BJNRZll [BJNRZ2] for further discussion of this and related 
topics. A further interesting direction for future study would be to obtain pointwise Green 
function bounds as in |JZN1 IMaZ3] for the metastable case, generalizing the bounds ob- 
tained in [OZl] for perturbations of unstable constant solutions. A related problem is to 
obtain rigorous nonlinear instability results in the metastable case. (Recall that essential 
instabilitiy does not immediately imply nonlinear instability due to the absence of a spectral 
gap.) 



A High Frequency Bounds 

Here, we record the quantities necessary to compute the high frequency bounds of Corol- 



lary 4.3 in terms of the original homoclinic orbit and of quantity a defined in (3.3) 
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